clear
A=[0 0 0 0;-1 0 0 0;0 0 0 1;0 -11.825 11.825 0];

B=[1.504 0;-0.471 0;0 0;0 0.671];
C=[-0.7 1 0 0;0 -1.207 1.207 0];
D=[0 0;0 -0.0165];
P=[-4.0+i;-4-i];
E=[1 0;0 1];
Q1=[-0.7 1 0 0;0 -1.207 1.207 0;1 0 1 0;0 1 0 1];
Q=inv(Q1);
A1=Q1*A*Q;
B1=Q1*B;
C1=C*Q;
q1=[Q(1:4,1:2)];q2=[Q(1:4,3:4)];
A11=[A1(1:2,1:2)];A12=[A1(1:2,3:4)];
A21=[A1(3:4,1:2)];A22=[A1(3:4,3:4)];
B11=[B1(1:2,1:2)];
B22=[B1(3:4,1:2)]; 
L=place(A22,A12,P)
Z1=A22-L*A12
Z2=A21-L*A11
Z3=B22-L*B11
q1+q2*L
q2
q1

